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ABSTRACT 

The disc formation mechanism of B[e] supergiants is one of the puzzling phenomena in 
massive star evolution. Rapid stellar rotation seems to play an important role for the non- 
spherically symmetric mass-loss leading to a high-density disc or ring-like structure of neu- 
tral material around these massive and luminous objects. The radial density and temperature 
structure as well as the kinematics within this high-density material are, however, not well 
studied. Based on high-resolution optical spectra of a sample of B[e] supergiants in the Mag- 
ellanic Clouds we especially searched for tracers of the kinematics within their discs. Besides 
the well-known [Ol] lines, we discovered the [Call] AA7291, 7324 lines that can be used as 
a complementary set of disc tracers. We find that these lines originate from very high-density 
regions, located closer to the star than the [O l] A5577 line-forming region. The line profiles of 
both the [O l] and the [Ca ll] lines indicate that the discs or rings of high-density material are 
in Keplerian rotation. We estimate plausible ranges of disc inclination angles for the sample 
of B [e] supergiants and suggest that the star LHA 1 20-S 22 might have a spiral arm rather than 
a disc. 
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1 INTRODUCTION 

B[e] supergiants (B[e]SGs) are evolved massive stars with large 
. amounts of circumstellar material in the form of atomic and molec- 
' ular gas as well as dust. The presence of gas is indicated by the 
numerous emission lines from both forbidden and permitted tran- 
sitions of atoms in different ionization stages. Molecular emission 
has been reported for many B[e]SGs: TiO emission has been de- 
tected at optical wavelengths ( Zickgraf et al. 1989. Torres et al., 
in prepar ation) and CO emission has been observed in the near- 
infrared (J McGregor. Hillier & Hvla ndl Il988l: iMorris et all Il996l : 
iLiermann et al. 2010; Muratore et a l .120 101) . The presence of dust is 
obvious from the strong infrare d excess emission in the photomet- 
ric and spec troscopic data (e.g. JZickgraf et al.ll986l : lBonanos et al.l 
l2009ll201(J) . and the composition of this dust shows indication for 
both oxygen (crystalline silicates) and carbon (polycyclic aromatic 
hydrocarbons, PAHs) chemistry ( Kastner et al. ;2010,). 



* Based on observations collected with the ESO 2.2m-telescope in La Silla, 
Chile, under programme 076.D-0609(A) 

t E-mail: aret@aai.ee (AA); kraus@sunstel.asu.cas.cz (MK); fmura- 
tore@carina.fcaglp.unlp.edu.ar (MEM); borges@on.br (MBF) 



The co-existence of a normal B-supergiant radiation driven 
wind and a large amount of molecular and dusty material sug- 
gests that the geometry of the circumstellar material cannot be 
spherically S3'mmetric, and polarimetric studies of B[e]SGs in the 
Magellanic Clouds (e .g.,[Magalhaes 1992; Magalhaes et al. 200^; 
iMelgarej o et al.'2001) and of Galactic B[e] stars and B[e]SG can- 
didates JOudmaiier & Drew.1999.) supported a flat, disc-like geom- 
etry. Though studied for many decades, the structure and geometry, 
as well as the formation history of the B[e]SGs' discs still remain 
unknown. 



Izickgraf et al.l ( I1985I) suggested that the discs are formed 
by a slow, high-density equatorial outflow. Such an equa- 
torially outflowing disc could resul t from ra pid stellar rota- 
tion via the bi-stabili ty mechanism JLamers & Pauldrach 1991J; 
IPelupessv. Lamers & Vink 2000 ) or from the newly f ound slow- 
wind solutions JCurg ,2004 ; .Cure. Rial & Cidalel l2005i) that pro- 
duce a slow but high-density wind in the equatorial region. 
Rapid stellar rotation could also be the reason for recombi- 
nation o^theionized wind material in the equatorial plane 
(e.g., iKrauj 1 200a) , resulting in large amounts of neutral hy- 
drogen close to the stellar surface, as was recently found for 
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Table 1. Parameters of the observed B[e]SGs. 



Object 



Sp. Type E{B - V) 



Tcff 
(lO^K) 



L 



Orientation References 









SMC 








LHA115-S 18 


BO 


0.40 


25 


3-4.6 


it pole-on 


(1) 


LHA115-S65 


B2-3 


0.15-0.20 


17 
LMC 


5.0 


edge-on 


(2) 


LHA 120-S 12 


B0.5 


0.20-0.25 


23 


2.2 


intermediate 


(2) 


LHA 120-S 22 


B0-B0.5 


0.25-0.30 


23-26 


6-9.5 


edge-on 


(2) 


LHA 120-S 73 


B8 


0.12 


12 


2.9 


it pole-on 


(2) 


LHA 120-S 111 




0.28 


15.5 


11.2 




(4) 


LHA 120-S 127 


B0.5 


0.25 


22.5 


13 


pole-on 


(3) 


LHA 120-S 134 


BO 


0.20-0.25 


26 


7.9 


pole-on 


(2) 

















References: (1) Izick 
lMcGregoreTd] l ll988l) 



IZickgrafetal.! l ll989l) : (2) Izickgraf et al] i ll986h : (3) Izickgraf et all h985l) : (4) 



some B[e] stars and B[e]SGs i Kraus & Borges FemandesI 120051 : 
iKraus. Borges Femandes & de Araijidl2007ll201ol) ^ 

However, for some individual B[e]SGs recent detailed inves- 
tigations revealed that the simple outflowing-disc scenario does 
not hold. The set of [Ol] AA5577, 6300 lines trace physically dis- 
tinct regions within the inner parts of the massive disc or ring 
whi ch are so dense that hyd rogen is already predominantely neu- 
tral JKraus et alj2007ll2010h . On the other hand, information on the 
density, temperature, and kinematics within the molecular parts of 
the disc can be obtained from the stre ngth and structure of the CO 
first-overtone bands (e.g., Carr et al. 



1993 1: IChandleretal.il 19931; 



nrst-overtone t)ancls (e.g ., Carr et a l.MlWJ , ICtiandler et al.MlWjl: 
C^ I1995I: iNaiita et al.l I1996.: .Kraus et al.ll2000l: iLiermann etal] 



2OI0I : iMuratore et al.ll2010l) . The fluxes, temperatures, and kine 



matics obtained from these emissions hint towards the presence 
of a detached ring (or detached disc-like st ructure) rather than a 
disc structure formed from a steady out flow dLiermann et al.bOld : 
iMuratore et alll2010l: iKraus et alj|2010l) . 

To study the formation mechanisms of B[e]SGs' discs it is 
important to combine the information of density, temperature and 
especially kinematics obtained from different sources, like opti- 
cal emission lines and infrared molecular bands that trace different 
regions within the circumstellar material. The forbidden emission 
lines play a very important role in this study, because they are op- 
tically thin (e.g., Os terbrock 1989) and hence id eal tracers for den- 
sity and temperature (see e.g.. lKraus et al.l2005n . and their profiles 
contain the full kinematical information of t heir formation region 
jlgnace & Brimevenl2006l : iKraus et alJuOlOl) . However, most for- 
bidden emission lines are formed in the low-density regions. The 
chance to find forbidden lines originating from the high-density 
disc complementing the well-known set of the [Ol] lines is rather 
low. Also, high-quality optical and near-infrared observations for 
most of the B[e]SGs are still lacking. 

To overcome the difficulties in studying the structure and for- 
mation history of the B[e]SGs' discs we started an observational 
campaign focused on the search for typical density and velocity 
tracers in the spectra of these stars. Here we report on our results 
based on high-resolution optical spectra obtained for a total of eight 
B[e]SGs in the Magellanic Clouds. In these spectra we discovered 
besides the well-known [O I] lines a set of forbidden emission lines, 
the [Can] AA7291, 7324 lines, as a new, valuable disc tracer. 



2 OBSERVATION AND REDUCTION 

Our sample consists of two B[e]SGs from the Small Magellanic 
Cloud and six from the Large Magellanic Cloud (SMC and LMC, 
respectively). The objects are listed in Table[T] together with the 
stellar parameters and possible disc orientations as found in the lit- 
erature. 

The observations have been carried out on two consecutive 
nights, on 2005 December 10 and 11. We obtained high-resolution 
optical spectra using the Fiber-fed Extended Range Optical Spec- 
trograph (FEROS) attached to the 2.2-m telescope at the European 
Southern Observatory in La Silla (Chile). 

FEROS is a bench-mounted Echelle spectrograph with two 
fibres, each of them covering a circular sky area of 2" in diam- 
eter. The wavelength coverage extends from 3600 A to 9200 A 
and the spectral resolution is i? = 55 000 (in the region around 
6000 A). We made our own data reduction, since the pipeline pre- 
sented some problems du ring that mission. Based on this, we used 
available MIDAS routines dVerschueren et al.l 19971) , including some 
modifications related to the FEROS characteristics. Thus, our data 
are reduced differentially relative to a master stellar and flat-field 
frame, considering the usual steps for echelle data reduction with 
the proper location and merge of the spectral orders. The spectra 
are rebinned taking into account the natural resolution of the de- 
tector. The telluric and barycentric velocity corrections have also 
been performed. From the sky lines we measured the accuracy of 
the wavelength calibration, which is better than 0.5 kms~^. 

To identify and delete cosmic rays from the spectra, we took 
a minimum of two spectra per star. The spectral lines did not show 
any variability in the individual exposures, so we added them up 
to achieve a better signal-to-noise ratio (S/N). Exposure times to- 
gether with the final achieved signal-to-noise ratios are listed in Ta- 
ble 2. The stellar radial velocities obtained from the central wave- 
length of the [O I] A6300 line are also given in that table. 



3 RESULTS 

The FEROS spectra of the B[e]SGs display a large number of emis- 
sion lines from permitted and forbidden transitions. Emission lines 
from the latter are especially valuable, because forbidden lines are 
optically thin, and their profiles reflect the kinematics within their 
formation region. One interesting set of lines has been identified 
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Table 2. Exposure times (number of exposures is given in parentheses) and 
signal-to-noise ratios at 5500 A. The radial velocities Vr were derived using 
the [0 1] A 6300. 



Object 


Exposure (s) 


Total S/N 


K(kms-i) 


LHA115-S 18 


1800 (2) 
900 (1) 


50 


146 


LHA115-S65 


900 (2) 


55 


189 


LHA 120-S 12 


900 (2) 


60 


290 


LHA 120-S 22 


900 (2) 


70 


289 


LHA 120-S 73 


600(1) 
900 (1) 


100 


261 


LHA 120-S 111 


900 (2) 


80 


267 


LHA 120-S 127 


450 (2) 


50 


259 


LHA 120-S 134 


1800(1) 
900 (1) 


70 


267 



as the forbidden emission lines of neutral oxygen, [O I] AA5577, 
6300, 6364. While the two lines AA6300, 6364 arise from the same 
upper energy level and hence are formed in the same region within 
the disc, the [O I] A5577 line emerges from r egions of higher den- 
sity (for details see Kraus et al. 20071 . uOlOh , i.e., from distances 
very close to the star. The ratio [O I] A6300/A5577 serves as an 
ideal indicator for temperature, density, and ionization fraction of 
the disc material (Kraus et al. 2007, 2010). Since the emission of 
these two individual lines arises from physically distinct disc re- 
gions, i.e., regions with different velocities, their line fluxes and 
profiles deliver complementary information. This makes these two 
lines a valuable tool in the study of both the structure and kinemat- 
ics of the discs, as was recently shown by Kraus et al. (2010) for 
the B[e]SG LHA 1 15-S 65 in the SMC. 

Unfortunately, for many stars of our sample the [O I] A5577 
line is quite weak or even absent (see Figs[T]and|2j. In addition, 
the S/N ratio of the observed spectra in the wavelength region of 
this line is rather low (see Table|2). This hampers severely a proper 
determination of the physical parameters of their discs and empha- 
sizes the need for additional lines that are capable of tracing the 
disc properties at different distances from the star. 

Upon inspection of the FEROS spectra we found that all 
our B[e]SGs display strong line emission from the Call AA8498, 
8542, 8662 infrared triplet. Examples of the FEROS spectra 
around these lines are shown in Fig.|3]for three stars of our sam- 
ple with different Call line profiles and different strengths in 
the adjacent hydrogen lines of the Paschen series. This emis- 
sion is known to arise from high-density regions, and circum- 
stellar discs are especially favoured locations. Emission in the 
Call infrared triplet has been reported for a diversity of objects: 
T Tau ri stars (e.g., Hamann & Persson 1992a; Kwan & Fischeg 
l201lh . Herbig A e/Be stars (e.g., Hamann & Persson 1992b i), classi- 
cal B e stars (e.g., Polidan & Peters^l 976:. BrioL 1981. : .Jaschek et al.l 



1 19881 : lAndrillat. Jaschek &Jascheklll988lll990l). and oth er, non- 
supergiant B[e] stars dBorges Femandes et al .11200 11120091) . 

Furthermore, we discovered strong line emission in [Ca II] 
AA7291, 7324 in all our objects. So far, the appearance of these 
forbidden em ission lines ha s been reported for only a few othe r 
disc sources ( lHamamj|l994l : lHartigan. Edwards & PiersonI 120041) . 
Nevertheless, the detection of these forbidden lines is of particu- 
lar interest, because the energy levels from which the [Ca II] lines 
emerge are coupled to those from which the Ca II infrared triplet 
lines arise. The presence or absence of the [Ca II] lines could thus 



be a suitable tracer for specific temperature and density conditions 
within the disc or wind. 

In Figs[T] and [2] we plot the line profiles of the calcium lines 
together with those of the [O I] lines and of Ha. To compare widths 
and shapes of lines from different stars we plotted identical lines on 
the same radial velocity scale. Stellar radial velocities as listed in 
Table[2]have been subtracted. The close vicinity of the Ca II infrared 
triplet lines to adjacent Paschen lines could result in line-blending. 
The distances of the individual Paschen lines Pa(16) A8502, Pa(15) 
A8545 and Pa(13) A8665, from the Call lines are +157.7kms"\ 
+ 115.5 km s^^ and +99.7 km s^^, respectively. However, in the 
high-resolution FEROS spectra the shape of the Ca II lines is not 
distorted by the Paschen lines for most of our stars (see Figs[T] |2] 
and[3j. One exception is LHA 120-S 22, in which the Paschen lines 
are extremely strong. 



3.1 Line profile diversity 

The appearance of the optical spectra of the B[e]SGs has been de- 
scribed in detail in the literature (e.g., Izickgraf et al.lll985L Il986l 
1989). However, former observations were limited to wavelengths 
shorter than 7000 A. This is the reason why the strong emission in 
both the Ca II infrared triplet and the forbidden [Ca II] lines was 
not reported earlier Hence, we refrain from repeating a detailed 
spectral description, and instead we focus on these newly detected 
emission lines and discuss their shape with respect to the [O I] lines. 

The different line widths of the forbidden lines are expected 
to reflect the kinematics of physically distinct regions. We note an 
increase in line width starting with the [O I] AA6300, 6364 lines, 
through the [Ol] A5577 line, to the [Call] AA7291, 7324 fines for 
all stars but LHA 120-S 127. In addition, we notice a simultaneous 
change in line profile. If the [O I] AA6300, 6364 lines have a single- 
peaked, Gaussian shape, the [Call] AA7291, 7324 lines show slight 
indications for double-peaks. If the [O I] AA6300, 6364 lines show 
slightly double-peaked profiles, the [Call] AA7291, 7324 lines are 
wide and clearly doubly-peaked. The width and shape of the [O I] 
A5577 line is intermediate. However, a clear classification of its 
profile shape is difficult for most of our objects, because this line 
is usually weak and noisy. The profile shapes and the measured (if 
possible) peak-separations in velocities are listed in Table[3] 

The lines of the Call infrared triplet are always (much) 
broader than the forbidden emission lines, and their line profiles 
show a diversity ranging from double (Fig.[T) to multiple-peaked 
(Fig.|2]l- Where possible, the peak-separation in velocity has been 
measured and listed in Table|4] The double-peaked profiles of these 
lines, however, do not resemble those structures that refer to the 
formation within a rotating or equatorially outflowing medium as it 
seems to be the case for the forbidden lines. Instead, they probably 
stem from absoiption within a high-density medium. 

The emission in the Ha line is the most prominent feature 
in all our spectra. Its strength ranges from about 25 to about 250 
times the continuum flux. Its profile can be single-peaked, double- 
peaked, or of P Cygni type. Its wings extend for all objects to large 
velocities. These high velocities, however, do not reflect real kine- 
matics but are most probably caused by electron scattering in a 
high-density wind or disc. 



3.2 Disc structure and kinematics 

The obvious diversity we see in the line profiles indicates that the 
formation regions for different lines must be physically distinct. 
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Figure 1. Line profile variety of stars seen pole-on (S 18. S 127) and edge-on (S 65, S 22). For details see text. 
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Figure 3. Ca II infrared triplet in LHA 120-S 22, LHA 120-S 111 and LHA 120-S 127. The spectr-a of LHA 120-S 1 1 1 and LHA I20-S 127 were shifted upwards 
for abetter viewing (by 11 and 21, respectively). The gap around 8530 A is a short wavelength region that is not covered by the red orders of FEROS. 

Table 3. Shapes and peak separations (km s~^) of forbidden calcium and oxygen lines. 



Object 


[Call] 729 


A 


[Call] 7324 A 


[O I] 5577 A 




[0 1] 6300 A 


[0 1] 6364 A 




Shape 


Sep. 


Shape 


Sep. 


Shape 


Sep. 


Shape 


Sep. 


Shape 


Sep. 


LHA115-S 18 


double 


10 


double 


10 


single 


_ 


single 


_ 


single 


_ 


LHA115-S65 


double 


38 


double 


37 


double? 


^ 


double 


n. r. 


double 


8 


LHA 120-S 12 


double 


42 


double 


45 


no line 


- 


single 


- 


single? 


- 


LHA 120-S 22 


double 


30 


double 


30 


double 


20 


double 


13 


double 


14 


LHA 120-S 73 


double 


11 


double 


12 


single 


- 


double 


6 


double 


6 


LHA 120-S 111 


double 


36 


double 


36 


double? 


- 


double 


n. r. 


double 


n. r 


LHA 120-S 127 


single 


- 


single 


- 


single 


- 


single 


- 


single 


- 


LHA 120-S 134 


double 


45 


double 


46 


broad single? 


- 


double 


13 


double 


13 



n. r - not resolved. 



This allows us to draw some (though preliminary) conclusions 
about the possible disc kinematics and density structure of the sam- 
ple stars. 

As was recently shown bv lKraus et"ai] J2007ll201(ih . the [O I] 
A5577 line arises from regions of much higher density and hence 
closer to the star than the [O I] AA6300, 6364 lines. Thus, be- 
sides information on temperature and density this set of lines pro- 
vides details on the kinematics within the two line-forming regions. 
For the B[e]SG star LHA 115-S 65 which is viewed edge-on the 
[O I] A5577 line is broader than the other two [O I] lines and with 
stronger indication of a double-peaked profile. Simultaneous mod- 
elling of the [O I] line profiles and line luminosities led to the con- 
clusion that the disc around LH A 115-S 65 is deta ched from the 
star and is in Keplerian rotation dKraus et alJuOlOh . Furthermore, 
this detached disc or ring has a rather small outflow velocity which 
decreases with distance from the star. These findings are in strict 
contrast to the general belief that B[e]SG stars' discs are the high- 
density eq uatorial outflow parts within the two-component wind 
scenario of Izickgraf et al.l(ll985h . 

Comparison of the [O I] line profiles for all our sample stars 
indicate that the A5577 line is indeed almost always the broadest 
one suggesting that Keplerian rotation rather than outflow might 



hold for all B[e]SG discs. Exceptions are the stars LHA 120-S 127 
that is viewed exactly pole-on so that all [Ol] lines are equally 
broad and LHA 120-S 12 for which the [O I] A5577 line could not 
be detected. If the discs around the B[e]SGs are indeed in Keplerian 
rotation, the broader and clearly double-peaked line profiles of the 
[Ca II] lines must originate from regions closer to the star than the 
[O I] A5577 line-forming region. Consequently, the density in the 
[Ca II] line-forming region must be higher. 

Recently, iHartigan et alj ^20041) investigated the [Call] 
A729 1/7324 line ratios in both the low- and the high-density limit 
and found that the ratio changes only a little from 1 .495 in the for- 
mer to 1.535 in the latter. They suggested the critical electron den- 
sity to be on the order of ~ 5 x lO'^ cm"''. To obtain the [Call] 
A729 1/7324 flux ratio, we measured the equivalent widths of the 
lines. We supposed that the real flux ratio does not differ substan- 
tially from the measured ratio of the equivalent widths because the 
lines have a separation of only 30 A. To confirm this, we checked 
the influence of the different continuum strengths and the effect 
of dereddening due to interstellar extinction by applying Kurucz 
model atmospheres and using the extinction values from the litera- 
ture as given in Table[T] We found that for all our stars corrections to 
the measured equivalent width ratios are less than 1%. The final val- 
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Table 4. Shapes and peak sepai'ations (km s ^) of Ca II infrared triplet. 



Object 


Call 8498 A 


Call 8542 A 


Call 8662 A 




Shape 


Sep. 


Shape 


Sep. 


Shape 


Sep. 


LHA115-S 18 


double 


n.r. 


double 


23 


double 


19 


LHA115-S65 


double 


33 


double 


26 


double 


31 


LHA 120-S 12 


multi 


- 


triple 


- 


triple 


- 


LHA 120-S 22 


multi 


- 


double 


28 


double 


28 


LHA 120-S 73 


triple 


- 


triple 


- 


triple 


- 


LHA 120-S 111 


multi 


- 


triple 


- 


triple 


- 


LHA 120-S 127 


double 


16 


double 


17 


double 


25 


LHA 120-S 134 


multi 


- 


triple 


- 


triple 


- 



n. r. - not resolved. 

ues for the [Call] A7291/7324 line ratios range from L33 to 1.47. 
All these value s fall below th e value of 1.495 for the low-density 
limit found by iHartigan et al.l |2004). Thus, we cannot determine 
the electron density in the [Call] line forming region from the ob- 
served flux ratios. Instead, we may use the values found within the 
[Or] A5577 line-forming region as lower limits for the electron 
densities within the [Call] line-forming regions. However, so far 
proper density determinations within the [O I] line-forming regions 
have been perform ed for only one object, namely LHA 115-S 65 
JKraus et alj2010l) . 

The SMC star LHA 1 15-S 65 was repor ted to be viewed edge- 
on jZickgraf et al] Il986l : JKraus et alj 120101) . Its [Call] lines are 
much broader than the [O I] lines and clearly double-peaked. The 
electron density within the [Ol] A 5577 line-forming region was 
found to be ric ~ 9 X lO" cm"^ JKraus et alj|20ia) . The range 
in rotational velocities over which the [Call] lines seem to be 
formed spreads from about 52kms~^ as measured from the ex- 
tent of the line wings to about 19kms^^ as obtained from the 
peak-separation of their profiles. Current studies of the evolutionary 
link of B[e]SGs to other phases in massive star evolution indicate 
that they could be evolving back bluewards, which means that they 
are in a post-red su pergiant (post-RSG) or post-yellow hypergiant 
(post-YHG) phase JMuratore et al.ll2OI0l Muratore et al. in prepa- 
ration). Applying Keplerian rotation for the stellar mass obtained 
for such a scenario and given in Table|5] we find that the velocity 
of 52 km s^^ corresponds to a distance of about 22 i?, and the ve- 
locity of I9kms^^ to about 160 7?,. Thus, most of the emission 
originates from the disc regions at roughly 160 i?*. If we assume 
that the disc has formed from stellar mass-loss, then its radial den- 
sity structure has a r~^ distribution according to the law of mass 
conservation. Consequently, the density increases by about a factor 
of > 50 over the [Ca II] emission region (from the outer to the inner 
edge) and becomes (much) higher than the critical density, meaning 
that we are in the high-density regime of the [Ca II] line formation. 

The rather large range in density might be one reason why the 
values of the [Ca II] flux ratios do not agree wi th either of the limit- 
ing values computed by IHartigan et alj ( 120041) . However, it should 
be noted that their values have been determined under the assump- 
tion that the material is ionized and free electrons are the dominant 
collision partners. Such a scenario does not hold for the discs of the 
B[e]SGs. Instead, the hydrogen ionization fractions were found to 
be very small within the [O I] line-forming regions (but still deliver- 
ing an electron density higher than the critical density for [Call]). 
The discs around B[e]SGs can t hus be c onsidered as being pre- 
dominantely neutral in hydrogen JKraus et al. 2007, 2010). Hence, 
collisions with hydrogen atoms certainly dominate over collisions 
with free electrons and need to be considered together with a proper 



radiation transfer in the lines of the infrared triplet. A proper mod- 
elling of the Ca II atom is beyond the scope of the current work. 
Nevertheless, we may claim that the [Ca II] lines provide a new im- 
portant and appropriate tracer for the structure and kinematics of 
B[e]SG stars' discs. 



3.3 Masses of B[e]SGs - their evolutionary phase 

The mass of LHA 115-S 65 used for the determination of the dis- 
tance of and the density within the [Call] line-forming region 
was obtained considering that the star is evolving bluewards in 
the Hertzsprung-Russell diagram (HRD). In such a scenario, it 
could have lost the mass that is now forming its disc during the 
passage through a RSG and/or YHG phase (see Fig.|4j. So far, 
not much is known about the real evolutionary phase of B[e]SGs 
and several suggestions about their present stage can be found in 
the literature. For instance, it has been s uggested that B[e]SGs 
are the results of massive binary mergers (Langer & Heger 1993; 
IPodsiadlowski. Morris & Iva nova 2006). As of now such a scenario 
seems to hold for only one object, the SMC stars R4. However, one 
should avoid to associate all B[e]SGs with stellar mergers, because 
the spectral characteristics of R4, its position in the HRD just at 
the end of the main sequence, and its strong photometric variabil- 
ity, which is not common for B[e]SGs, displaces R4 from the rest 
of the B[e]SG sample. We therefore consider the merger scenario 
as the least likely common origin for the B[e] phenomenon in the 
supergiant stars. 

Some B[e]SGs were suggested to be evolutionary linked to 
Luminous Blue Variables (LBVs) with which they share about the 
same region in the HRD (see Fig.|4j. For instance, t he star V39 in 
the low metallicity galaxy IC 1613 was classified by Herrero et al.l 
(2010) as either a B[e]SG or an L BV candidate, and the B[e]SG star 
LHA 1 15-S 65 was classified bv iKraus et alj J2010l) as a pre-LBV. 
The LBVs are luminous objects (log L/Lq Js 5.4). They were 
found to separate into two groups, and this separation happens at 
log I//I/Q ~ 5.8. While the more luminous LBVs were suggested 
to be post-main sequence stars evolving redwards, their lower lu- 
minosity counterparts are thought to evolv e bluewards, i.e., they are 
sugge sted to be post-RSG s or post- YHGs jHumphrevs & DavidsonI 
1994: lMevnet et alj201 ll. see also Fig.gJ. 

Half of our objects certainly fall into this low-luminosity 
regime, while two others are just at the border, but could still fit to 
the same scenario because their location in the HRD indicates that 
they might have evolved from a 40-42 Mq progenitor star, which 
is about the limiting m ass between the two evolutionary scenarios 
for LBVs described bv lMevnet et al.l ( 1201 ih . The possible classifi- 
cation of at least one of our sample B[e]SGs as a pre-LBV further 
suggests that the bluewards evolution in this low-luminosity regime 
could indeed go from the Y HG over th e B[e]SG to the LBV phase. 

In recent investigations |Kastneret al. ( 2010) studied the dusty 
discs of the Magellanic Cloud B[e]SGs. They found convincing 
analogies to the discs around lower mass stars in a post-asymptotic 
giant branch phase, so they concluded that the B[e]SGs are most 
likely the more massive counterparts in a post-RSG (or post-YHG) 
phase. And further support for such a scenario comes from re- 
cent investigations of the ^"^C footprint for the two least massive 
B[e]SGs in our sample (LHA 120-S 73 and LHA 120-S 12) reveal- 
ing that a post-YHG scenario is most likely for LHA 120-S 73 
( lMuratoreetaLll2 010'). The ^^C/^C ratio found for LHA 120-S 12 
was in agreement only with an evolved nature of the star. Whether 
LHA 120-S 12 could also be in a post-YHG phase needs to be tested 
based on predictions of the ^^C/^'^C ratio from evolutionary mod- 
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Fig ure 4. Our B [e]SG sample in comparison to other massive post-main sequenc e phases. Stellar param eters are taken from lHumphrevs & Davidsonl 41994) 
and lGroh et"aD f2009) for the LBVs, from de Jaaer (1998) for the YHGs, and from lLevesgue et alj I2003) for the RSGs. Stellar evolutionary tracks for rotating 
stars at LMC metallicity from lMevnet & Maeder t2005.) are overplotted. 



els for rotating stars at LMC metallicity extend ing to lower initial 
masses than provided by the current models of iMevnet & Maedetl 
feOOS) . but such models are not yet available. 

However, a post- YHG scenario certainly fails for the most lu- 
minous B[e]SGs, for which no YHG phase exists. At which stage 
in the evolution they might be (just leaving the main-sequence as 
was suggested for the most luminous LBVs or on a blue loop) is 
not known. 



3.4 Disc inclination angles 

We now return to the [Ca II] lines. The measured peak separation 
of their line profiles (see Table|3) can be used to estimate possible 
ranges in disc inclination angles. For this, we need to know the ap- 
propriate distance from the stars at which the bulk of the [Ca II] 
lines form. For LHA 115-S 65 we found that this distance is about 
160 Rf The density reached at that distance should be very similar 
for the [Ca II] line-forming region for all B[e]SGs. We can estimate 
these distances based on the disc mass fluxes that were determined 
by[zickgraf (1998) for all our sample stars but LHA 120-S 111. 
Both the ratios of disc mass fluxes of the sample stars with respect 
to LHA 115-S 65 and the resulting distances at which the required 
density is reached are given in Table|5] 

With the obtained distances, inclination angles can be esti- 
mated from the measured peak-separations together with the as- 
sumption of a Keplerian rotation law. To calculate Keplerian rota- 
tion velocities the current stellar masses need to be known. In Fig.|4] 



we show the positions of our objects in th e HRD together with evo- 
lutionary tracks from Mev net & Maedeii (|2005') computed for ro- 
tating stars at LMC metallicity. The evolutionary tracks for SMC 
metallicity do not differ substantially from those for the LMC, so 
we refrain from plotting both sets of tracks. As discussed in the 
previous section, the B[e]SGs with \ogL/ Lq < 5.8 might be in a 
post- YHG phase, while for the most luminous ones it is not clear 
whether they are evolving to the red or to the blue side of the HRD. 

Assuming that B[e]SGs are evolving bluewards and interpo- 
lating available evolutionary tracks, we obtained zero-age main- 
sequence and current masses for all objects as listed in Table[5] For 
the most massive stars in the sample we also list current masses ob- 
tained under the assumption that these stars might just have evolved 
off the main-sequence. The stellar radii have been obtained from 
the temperature and luminosity values given in Table[T] 

The line-of-sight velocities, t^ios, result from the peak- 
separations of the [Ca II] lines. If all stars were seen edge-on, these 
velocities would correspond to the Keplerian velocities. To get an 
idea about maximum distances of the [Ca II] line-forming regions 
from the central stars we calculated the radial distances resulting 
from the peak-separation velocities. These values are listed in Ta- 
ble|5]in terms of stellar radii. Obviously, these distances are quite 
large for those stars that are considered to be seen more or less 
pole-on like, e.g., LHA 115-S 18. 

We computed the Keplerian rotation velocities for all objects 
at the distance r([Ca II]) where the bulk of the [Ca II] line emis- 
sion is thought to originate. From the ratio of the observed rota- 
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Table 5. Stellar m asses, Keplerian rotation velocities, and inclination angles as estimated from the [Ca II] line profiles. Values of the disc mass fluxes, 
fa, are taken from Izickgrj l l 19981) . 



Object AfzAMS Mcurront Rt /d//d(S65) r-([Ca II]) VKcpir) vias r(vKap = I'los) suggested i 

(Mq) (Mq) {Rq) (R*) (kms-i) (kms-i) (/?,) (°) 



LHA115-S 18 


29 ±2 


19 ±3 


33 ±3 


2.15 


235 


21.7 ±2.7 


5 


4491 ± 1104 


13.5 ±2 


LHA115-S65 


35 ±1 


25 ±1 


82 


1.00 


160 


19.0 ±0.4 


19 


161 ±7 


84 ±6 


LHA 120-S 12 


22 ±2 


14 ±2 


30 


2.15 


235 


19.8 ±1.0 


21-22.5 


191 ±41 


~90 


LHA 120-S 22 


41 ±4 


25 ±4 


49 ±11 


3.85 


314 


18.0 ±3.5 


15 


469 ± 173 


67 ±23 




49 ±4=' 


41 ±4 








23.1 ±3.7 




763 ± 240 


42 ±8 


LHA 120-S 73 


25 ±1 


15 ±1 


125 


0.78 


141 


12.7 ±0.4 


6 


634 ± 42 


28 ±1 


LHA 120-S 1 1 1 


60 ±5 
60 ±5=^ 


40 ±5 
48 ±5 


147 


unknown 


- 


- 


18 


171 ± 12 
192 ± 20 




LHA 120-S 127 


^60 


unknown 


75 


1.31 


183 


- 


unresolved 


- 


~0 


LHA 120-S 134 


42 ±1 


26 ±1 


44 


2.46 


251 


21.2 ±0.5 


23 


196 ±8 


~90 




50 ±2=" 


42 ±2 








26.9 ±0.6 




344 ± 16 


59 ±2 



' assuming the star has just left the main-sequence. 



tional velocity projected to the line-of-sight and the expected ro- 
tational velocity, we obtained a range in inclination angles. These 
values are included in Table|5] Good agreement with formerly sug- 
gested orientations is achieved for most objects. In particular, the 
stars LHA 1 15-S 18 and LHA 120-S 73 are found to be ±pole-on. 

The inclination angle of LHA 120-S 111 has not yet been de- 
termined. The position of this star in the HRD shows that it belongs 
to the most luminous objects. The masses estimated for both evo- 
lutionary scenarios can only be considered as rough estimates, and 
the missing value of its disc mass flux hampers a proper inclination 
determination. Nevertheless, considering the mass ranges as listed 
in Table|5] maximum distances of the [Ca II] line-forming region as 
obtained for a possible edge-on orientation are only slightly larger 
than the one obtained for LHA 1 15-S 65. The line strengths of the 
[Ol] and [Call] lines arecomparable to those of the other B[e]SGs, 
so that its disc mass flux should not be considerably smaller com- 
pared to LHA 1 15-S 65. Hence, we might conclude that LHA 120- 
S 1 1 1 could be oriented close to edge-on. 

The values obtained for the B[e]SGs at the border between 
low- and high-luminosity regime, LHA 120-S 22 and LHA 120- 
S 134, differ significantly for the two different evolutionary scenar- 
ios; while a close to edge-on orientation is obtained for the post- 
YHG scenario, an intermediate inclination angle is obtained for the 
redwards evolutionary scenario. 

For LHA 120-S 134, both results clearly dis agree with the 
pole-o n orientation (see Table[T) suggested by Izickgraf et al.l 
( 11986 1). The broad and clearly double-peaked line profiles of its 
[Call] lines definitely speak against a pole-on orientation (see 
Fig.|2j. Which scenario will turn out to be the correct one needs 
further detailed investigation. 

The orientation of LHA 120-S 22 in the redwards evolutionary 
scenario delivered a value of only i 2^ 42 ± 8°. But t his object has 



been reported to be viewed edg e-on ( Zickgraf et al.ll986) . This ori- 
entation has been confirmed bv lChu et al.l ( 1200 3h who discovered a 
reflection nebula in the vicinity of LHA 120-S 22 that mirrors dif- 
ferent viewing angles. Interestingly, a (close-to) edge-on orienta- 
tion is found in the post-YHG scenario. 

Inspection of the line profiles of LHA 120-S 22 reveals that 
this object has another peculiarity: both sets of forbidden lines are 
double-peaked and clearly asymmetric. But while the [Ca II] lines 
show a much stronger red peak, the [O I] lines display a stronger 
blue peak. Considering that the [Ca II] lines are supposed to orig- 
inate from distances closer to the star than the [O I] lines, this be- 



haviour indicates the presence of inhomogeneities in the circum- 
stellar material like a density wave in the disc or a spiral arm-like 
structure. Such structures have recently been seen on images of the 
disc of a young stellar objecij. So, LHA 120-S 22 might be the first 
B[e]SG star with spiral arms. In any case, it is worth to study both 
the kinematics of its circumstellar material and the time variability 
of the line profiles in more detail. 

We would like to emphasize that the inclination angles listed 
in Table[5]can only be considered as rough estimates, because sev- 
eral assumptions have been made for their derivation. The one with 
the highest uncertainty is the distance at which the [Ca II] line emis- 
sion forms, because we implicitly adopt that the ionization fraction 
found for LHA 1 15-S 65 is the same in all B[e]SG star discs. This 
must not be true. A higher value for the ionization fraction implies 
a higher electron density, and the distance at which the [Ca II] emis- 
sion forms is shifted farther out. The Keplerian rotation velocities 
are lower there and consequently the inclination of the disc would 
be larger. Proper values for the ionization fractions in the discs of 
all B[e]SGs are therefore needed. This emphasizes the need of flux 
calibrated data for all the sample stars. 



4 CONCLUSIONS 

We report on the discovery of the [Call] AA7291, 7324 lines in 
high-resolution spectra of a sample of B[e] supergiants. These lines 
originate from different regions than the already well-known disc 
tracers, the [Ol] lines. The [Call] lines trace regions of higher 
density and thus must form closer to the star. Combination of the 
kinematical information obtained from the line profiles of the [O I] 
and [Call] allows to suggest that the material around the B[e] 
supergiants is in Keplerian rotation. We estimated disc inclina- 
tions which are in reasonably good agreement with former guesses. 
However, the clearly asymmetric line profiles seen in the spectra of 
the star LHA 120-S 22 indicate that the material is not in a disc or 
ring-like structure. We suggest this star might have a spiral arm 
seen edge-on. 

In addition, we discuss plausible evolutionary phases for the 
B[e]SGs. Although perhaps for some individual B[e]SGs (like R4 
in the SIVIC) a different scenario might hold, a blueward evolution 
of (most of) the B[e]SGs seems to be likely. As suggested bv lKrauj 



http://www.nasa.gov/topics/universe/features/possible-planets.html 
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2009 ) and proven by iLiermann et al.l (120101) and iMuratore et alj 



20101), the C footprint is a powerful tool to determine the evolu- 



tionary phase of B[e]SGs. Hence reliable data of the C footprint 
for all B[e]SGs are needed, and we recently started an observational 
campaign to study the ^^C enrichment in the discs of all B[e]SGs. 
Results from these observations will certainly provide a better ba- 
sis for a more detailed discussion about the evolutionary phase of 
B[e]SGs and their possible links to other evolved phases of massive 
star evolution. 
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